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We present a brief review of various approaches to late time acceleration of universe. The cosmo- 
logical relevance of scaling solutions is emphasized in case of scalar held models of dark energy. The 
underlying features of a variety of scalar held models is highlighted. Various alternatives to dark 
energy are discussed including the string curvature corrections to Einstein-Hilbert action, higher 
dimensional effects, non-locally corrected gravity and f(R) theories of gravity. The recent develop- 
ments related to f(R) models with disappearing cosmological constant are reviewed. 
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I. INTRODUCTION 

Os . 
O ■ 

The accelerated expansion has played a very important role in the history of our universe. Universe is believed to 
£Nl ' have passed through inflationary phase at early epochs and there is a growing faith that it is accelerating at present. 
The late time acceleration of the universe, which is directly supported by supernovae observations, and indirectly, 
through observations of the microwave background, large scale structure and its dynamics, weak lensing and baryon 
, oscillations, poses one of the most important challenges to modern cosmology. 

■ Einstein equations in their original form, with an energy-momentum tensor for standard matter on the right hand 
side, cannot account for the observed accelerated expansion of universe. The standard lore aimed at capturing this 
important effect is related to the introduction of the energy-momentum tensor of an exotic matter with large negative 
pressure dubbed dark energy in the Einstein equations. The simplest known example of dark energy (for recent 
reviews, see [l|) is provided by the cosmological constant A. It does not require adhoc assumption for its introduction, 

'■ as is automatically present in the Einstein equations, by virtue of the Bianchi identities. 

The field theoretic understanding of A is far from being satisfactory. Efforts have recently been made to obtain 
l— ~~ 1 ■ A in the framework of string theory, what leads to a complicated landscape of de-Sitter vacua. It is hard to believe 
I | that we happen to live in one of the 10 500 vacua predicted by the theory. One might take the simplified view that, 

■ like G, the cosmological constant A is a fundamental constant of the classical general theory of relativity and that it 
\Q ' should be determined from large scale observations. It is interesting to remark that the ACDM model is consistent 

. with observations at present. Unfortunately, the non-evolving nature of A and its small numerical value lead to a 
non-acceptable fine-tuning problem. We do not know how the present scale of the cosmological constant is related to 
Planck's or the supersymmetry breaking scale; perhaps, some deep physics is at play here that escapes our present 
t-H ■ understanding. 

The fine-tuning problem, associated with A, can be alleviated in scalar field models which do not disturb the 
0\ . thermal history of the universe and can successfully mimic A at late times. A variety of scalar fields have been 
investigated to this end[l[; some of them are motivated by field/string theory and the others are introduced owing to 
. phenomenological considerations. It is quite disappointing that a scalar field description lacks predictive power; given 
a priori a cosmic evolution, one can always construct a field potential that would give rise to it. These models should, 
however, not be written off, and should be judged by the generic features which might arise from them. For instance, 
the tracker models have remarkable features allowing them to alleviate the fine-tuning and coincidence problems. 
Present data are insufficient in order to conclude whether or not the dark energy has dynamics; thus, the quest for 
the metamorphosis of dark energy continues Q 

One can question the standard lore on fundamental grounds. We know that gravity is modified at small distance 
scales; it is quite possible that it is modified at large scales too where it has never been confronted with observations 
directly. It is therefore perfectly legitimate to investigate the possibility of late time acceleration due to modification of 
Einstein-Hilbert action. It is tempting to study the string curvature corrections to Einstein gravity amongst which the 
Gauss-Bonnet correction enjoys special status. A large number of papers are devoted to the cosmological implications 
of string curvature corrected gravity 0, 0, H, 0, S S S 03 EH Q3- These models, however, suffer from several 
problems. Most of these models do not include tracker like solution and those which do are heavily constrained by the 
thermal history of universe. For instance, the Gauss-Bonnet gravity with dynamical dilaton might cause transition 
from matter scaling regime to late time acceleration allowing to alleviate the fine tuning and coincidence problems. 
But it is difficult to reconcile this model with nucleosynthesis [H, Q constraint. The large scale modification may also 
arise in extra dimensional theories like DGP model which contains self accelerating brane. Apart from the theoretical 
problems, this model is heavily constrained by observation. 

On purely phenomenological grounds, one could seek a modification of Einstein gravity by replacing the Ricci scalar 
by generic function f (j?)(l3lll4l [l5|. The f(R) gravity theories giving rise to cosmological constant in low curvature 
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regime are faced with difficulties which can be circumvented in f(R) gravity models proposed by Hu-Sawicki and 
Starobinsky [l6|, [f7| (see Ref.[lE| on the similar theme). These models can evade solar physics constraints by invoking 
the chameleon mechanism |l6t If 7l. Il9j . An important observation has recently been made in Refs.[2(], [2l[(see also 
Ref. [12] on the related theme), namely, the minimum of scalaron potential which corresponds to dark energy can be 
very near to <f> — or equivalently R — oo. As pointed out in Ref.[19j], the minimum should be near the origin for 
solar constraints to be evaded. Hence, it is most likely that we hit the singularity if the parameters are not properly 
fine tuned. This may have serious implications for relativistic stars [23j| - 
In what follows we shall briefly review the aforesaid developments. 



II. LATE TIME ACCELERATION AND COSMOLOGICAL CONSTANT 



Einstein equations exhibit simple analytical solutions in a homogeneous and isotropic universe. The dynamics in 
this case is described by a single function of time a(t) dubbed scale factor, 

ff2 _ a 2 _ %ixGp K 
a 2 3 a 2 
p + 3H(p + p) = 0, 

where p is designates the total energy density in the universe. Three different possibilities, K — 0, K > or K < 
correspond to flat geometry, hyperbolic geometry and geometry of the sphere correspondingly. Evolution can not 
change the nature of a particular geometry. What geometry we live in, depends upon the energy content of the 
universe, 

*=H 2 (fl(f)-l) 

n = p/p c , p c = 3H 2 /8irG 

Observations on CMB reveal that we live in a nearly critical universe, K = or p = p c which is consistent with 
inflationary paradigm. The equation for acceleration has the following form, 

a 4ttG 

- = — 5- (P + m 

a 3 

a > p < — — : DarkEnergy. 
3 

Thus an exotic fluid with large negative pressure is needed to fuel the accelerated expansion of universe. Let us note 
that pressure corrects the energy density and positive pressure adds to deceleration where as the negative pressure 
contributes towards acceleration. It might look completely opposite to our intuition that highly compressed substance 
explodes out with tremendous impact whereas in our case pressure acts in the opposite direction. It is important 
to understand that our day today intuition with pressure is related to pressure force or pressure gradient. In a 
homogeneous universe pressure gradients can not exist. Pressure is a relativistic effect and can only be understood 
within the frame work of general theory of relativity. Pressure gradient might appear in Newtonian frame work in 
an inhomogeneous universe but pressure in FRW background can only be induced by relativistic effects. Strictly 
speaking, it should not appear in Newtonian gravity; its contribution is negligible in the non-relativistic limit. Indeed, 
in Newtonian cosmology, acceleration of a particle on the surface of a homogeneous sphere with density p and radius 
R is given by, 

R = -—G P R (1) 

The simplest possibility of dark energy is provided by cosmological constant which does not require an adhoc assump- 
tion for its introduction; it is automatically present in Einstein equations by virtue of Bian'chi identities, 

G M „ = R^v - -j9iivR = SttGT^ - Ag^ (2) 



The evolution equations in this case become, 



" 2 = ^+f < 3 > 

a 4irG . A 
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In case the universe is dominated by A, it follows from the continuity equation that pa — ~PA and Eq.Q tells us 
that a positive cosmological constant contributes to acceleration. 
Observations of complimentary nature reveal that, 

• n tot ~ 1 - CMB, 

• fi m ~ 0.3 -Large scale structure and its dynamics, 

• Ude — 0.7 - high redshift la Supernove, 

which is independently supported by data on baryon oscillations and weak lensing. 

Observations at present do not rule out the phantom dark energy with w < — 1 corresponding to super acceleration. 
In this case the expanding solution takes the form, 



a(t) = (t s -t) n , (n = 2/3(1 + w)) 



(5) 



where t s is an integration constant. It is easy to see that phantom dominated universe will end itself in a singularity 
known as big rip or cosmic doomsday, 



H 



t,-t 



R = 6 



6 



n(n — 1) + n 2 



(6) 
(7) 



The big rip singularity is characterized by the divergence of H and consequently of the curvature after a finite time in 
future. It should be noted that when curvature becomes large, we should incorporate the higher curvature corrections 
to Einstein- rlilbert action which can modify the structure of the singularity [13] • Big rip can also be avoided in specific 
models of phantom energy [25j|. 



A. Issues associated with A 



There are important theoretical issues related to cosmological constant. Cosmological constant can be associated 
with vacuum fluctuations in the quantum field theoretic context. Though the arguments are still at the level of 
numerology but may have far reaching consequences. Unlike the classical theory the cosmological constant in this 
scheme is no longer a free parameter of the theory. Broadly the line of thinking takes the following route. 

The quantum effects in GR become important when the Einstein Hilbert action becomes of the order of Planck's 
constant; this happens at the Planck's length Lp — 10~ 32 cm corresponding to Planck energy which is of the order 
of ~ 10 72 GeU 4 . The ground state energy dubbed zero point energy or vacuum energy of a free quantum field 
is infinite. This contribution is related the ordering ambiguity of fields in the classical Lagrangian and disappears 
when normal ordering is adopted. Since this procedure of throwing out the vacuum energy is adhoc, one might 
try to cancel it by introducing the counter terms. The later, however requires fine tuning and may be regarded as 
unsatisfactory. The divergence is related to the modes of very small wave length. As we are ignorant of physics around 
the Planck scale, we might be tempted to introduce a cut off at L p and associate with this a fundamental scale. Thus 

we arrive at an estimate of vacuum energy py ~ Af 4 (corresponding mass scale- My ~ p]j 4 ) which is away by 120 
orders of magnitudes from the observed value of this quantity. The vacuum energy may not be felt in the laboratory 
but plays important role in GR through its contribution to the energy momentum tensor as < T^„ >o= —pv9nv 
(py = A/8irG, Mp — l/87rG)and appears on the right hand side of Einstein equations. 

The problem of zero point energy is naturally resolved by invoking supersymmetry which has many other remarkable 
features. In the supersymmetric description, every bosonic degree of freedom has its Fermi counter part which 
contributes zero point energy with opposite sign compared to the bosonic degree of freedom thereby doing away with 
the vacuum energy. It is in this sense the supersymmetric theories do not admit a non-zero cosmological constant. 
However, we know that we do not leave in supersymmetric vacuum state and hence it should be broken. For a viable 
supersymmetric scenario, for instance if it is to be relevant to hierarchy problem, the suppersymmetry breaking scale 
should be around M susy ~ 10 3 GeV. We are still away from the observed value by many orders of magnitudes. At 
present we do not know how Planck scale or SUSY breaking scales is related to the observed vacuum scale. 

At present there is no satisfactory solution to cosmological constant problem. One might assume that there is 
some way to cancel the vacuum energy. One can then treat A as a free parameter of classical gravity similar to 
Newton constant G. However, the small value of cosmological constant leads to several puzzles including the fine 
tuning and coincidence problems. The energy density in radiation at the Planck scale is of the order of 10 72 GeV 4 
or pa/ Pr ~ 10~ 120 Thus the vacuum energy density needs to be fine tuned at the level of one part in 10~ 120 around 
the Plank epoch, in order to match the current universe. Such an extreme fine tuning is absolutely unacceptable 
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FIG. 1: Desired evolution of field energy density p$ (pB is the background energy density). The field energy density in case 
of undershoot and overshoot joins the scaling solution for different initial conditions. At late times, the scalar field exits the 
scaling regime to become the dominant component. 




FIG. 2: Evolution of p<£ and ps in absence of the scaling solution. The scalar field after its energy density overshoots the 
background gets into locking regime and waits till its energy density becomes comparable to p_g. It then begin to evolve and 
over takes the background, similar picture holds in case of the overshoot. 



at theoretical grounds. Secondly, the energy density in cosmological constant is of the same order of matter energy 
density at the present epoch. The question what causes this coincidence has no satisfactory answer. 

Efforts have recently been made to understand A within the frame work of string theory using flux compactification. 
String theory predicts a very complicated landscape of about 10 500 de-Sitter vacua. Using Anthropic principal, we 
are led to believe that we live in one of these vacua! It is easier to believe in God than in these vacua! 



III. SCALAR FIELD DYNAMICS RELEVANT TO COSMOLOGY 



The fine tuning problem associate with cosmological constant led to the investigation of cosmological dynamics of a 
variety of scalar field systems such as quientessence, phantoms, tachyons and K-essence. In past years, the underlying 
dynamics of these systems has been studied in great detail. It is worthwhile to bring out the broad features that makes 
a particular scalar field system viable to cosmology. The scalar field model aiming to describe dark energy should 
possess important properties allowing it to alleviate the fine tuning and coincidence problems without interfering with 
the thermal history of universe. The nucleosynthesis puts an stringent constraint on any relativistic degree of freedom 
over and above that of the standard model of particle physics. Thus a scalar field has to satisfy several important 
constraints if it is to be relevant to cosmology. Let us now spell out some of these features in detail. In case the 
scalar field energy density p^ dominates the background (radiation/matter) energy p#, the former should redshift 
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faster than the later allowing radiation domination to commence which in tern requires a steep potential. In this case, 
the field energy density overshoots the background and becomes subdominant to it. This leads to the locking regime 
for the scalar field which unlocks the moment the p^ is comparable to ps- The further course of evolution crucially 
depends upon the form the potential assumes at late times. For the non-interference with thermal history we require 
that the scalar field remains unimportant during radiation and matter dominated eras and emerges out from the 
hiding at late times to account for late time acceleration. To address the issues related to fine tuning, it is important 
to investigate the cosmological scenarios in which the energy density of the scalar field mimics the background energy 
density. The cosmological solution which satisfy this condition are known as scaling solutions, 

— = const. (8) 

PB 

The steep exponential potential V(4>) ~ exp(X(f>/Mp) with A 2 > 3(1 + wb) in the frame work of standard GR gives rise 
to scaling solutions. Nucleosynthesis further constraints A. The introduction of a new relativistic degree of freedom at 
a given temperature changes the Hubble rate which crucially effects the neutron to proton for temperature of the order 
of one MeV when weak interactions freeze out. This results into a bound on A, namely, Sl^ = 3(l + u>s)/A 2 < 0.13 — 0.2 
or A > 4.5. In this case, for generic initial conditions, the field ultimately enters into the scaling regime, the attractor 
of the dynamics and this allows to alleviate the fine tuning problem to a considerable extent. The same holds for the 
case of undershoot, see FigQ] 

Scaling solutions, however, are not accelerating as they mimic the background (radiation/matter). One therefore 
needs some late time feature in the potential. There are several ways of achieving this: (1) The potential that mimics 
a steep exponential at early epochs and reduces to power law type V ~ 2p at late times gives rise to accelerated 
expansion for p < 1/2 as the average equation of state < >= (p — l)/(p+l) < —1/3 in this case 26]. (ii) The 
steep inverse power law type of potential which become shallow at large values of the field can support late time 
acceleration and can mimic the background at early times [27|. 

The solutions which exhibit the aforesaid features are referred to as tracker solutions. For a viable cosmic evolution 
we need a tracker like solution. 

Recently, a variety of scalar field models such as tachyon and phantom were investigated as candidates of dark 
energy. In case of tachyon with equation of state parameter ranging from —1 to 0, there exists no scaling solution 
which could mimic the realistic background (radiation/matter). Scaling solution which are possible in this case are 
associated with negative equation of state and are not of interest. In case of phantom scalar fields (scalar fields with 
negative kinetic energy) , there is no fixed point corresponding to scaling solution. These scenarios suffer from the fine 
tuning problem; dynamics in this case acquires dependence on initial conditions [2a] (see Figf5]). 

The second approach to late tim acceleration is related to the modification of left hand side of Einstein equations 
or the geometry of space time. In the past few years, several schemes of large scale modifications have been actively 
investigated. In what follows, we shall briefly describe the modified theories of gravity and their relevance to cosmology. 



IV. MODIFIED THEORIES OF GRAVITY AND LATE TIME ACCELERATION 



In view of the above discussion, it is perfectly legitimate to investigate the possibility of late time acceleration due 
to modification of Einstein-Hilbert action. Some of these modifications are inspired by fundamental theories of high 
energy physics where as the others are based upon phenomenological considerations. 



A. String curvature corrections 



It is interesting to investigate the string curvature corrections to Einstein gravity amongst which the Gauss-Bonnet 
correction enjoys special status. These models, however, suffer from several problems. Most of these models do not 
include tracker like solution and those which do are heavily constrained by the thermal history of universe. For 
instance, the Gauss-Bonnet gravity with dynamical dilaton might cause transition from matter scaling regime to late 
time acceleration allowing to alleviate the fine tuning and coincidence problems. Let us consider the low energy 
effective action, 



S = J <Fx^j 
- V(<j>) -f{<f)R\ 



1 



GB 



(9) 
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where Rq B is the Gauss-Bonnet term, 

R 2 GB = R 2 - \R ilv R^ + R a ^ v R a ^ v (10) 

The dilaton potential V((f>) and its coupling to curvature /(</>) are given by, 

V(<f>)~e^\ /(0)~ e -<"« (11) 

The cosmological dynamics of system (j9)) in FRW background was investigated in Ref.JH, It was shown that 
scaling solution ca be obtained in this case provided that // = A. In case fi =/= A, we have the de-Sitter solution. Hence, 
the string curvature corrections under consideration can give rise to late time transition from matter scaling regime. 
Unfortunately, it is difficult to reconcile this model with nucleosynthesis [f| [^constraint. 

B. DGP model 

In DGP model, gravity behaves as four dimensional at small distances but manifests its higher dimensional effects at 
large distances. The modified Friedmann equations on the brane lead to late time acceleration. The model has serious 
theoretical problems related to ghost modes superluminal fluctuations. The combined observations on background 
dynamics and large angle anisotropies reveal that the model performs worse than AC DM 29]. 

C. Non-local cosmology 

An interesting proposal on non-locally corrected gravity involving a function of the inverse d'Almbertian of the 
Ricci scalar, /(□ _1 i?)), was made in Refs. (30j For a generic function f(\3~ 1 R) ~ exp(aD~ 1 R) , the model can lead to 
de-Sitter solution at late times. The range of stability of the solution is given byl/3<a<2/3 corresponding to the 
effective EoS parameter w c fF ranging as oo < w e g < —2/3. For 1/3 < a < 1/2 and 1/2 < a < 2/3, the underlying 
system is shown to exhibit phantom and non-phantom behavior respectively; the de Sitter solution corresponds to 
a = 1/2. For a wide range of initial conditions, the system mimics dust like behavior before reaching the stable 
fixed point at late times. The late time phantom phase is achieved without involving negative kinetic energy fields. 
Unfortunately, the solution becomes unstable in presence of the background radiation/matter [30j. 

D. f(R) theories of gravity 

On purely phenomenological grounds, one could seek a modification of Einstein gravity by replacing the Ricci scalar 
by f(R)- The f(R) gravity theories giving rise to cosmological constant in low curvature regime are plagued with 
instabilities and on observational grounds they are not distinguished from cosmological constant. The recently intro- 
duced models of /(i?) gravity by Hu-Sawicki and Starobinsky (referred as HSS models hereafter) with disappearing 
cosmological constant [la. [l7j have given rise to new hopes for a viable cosmological model within the framework of 
modified gravity. The action of f(R) gravity is given bv[l3j . 



S ■■ 

which leads to the following modified equations 



(12) 



f'Ruv - V^f + (nf ~ i/) = SttGTV (13) 

which are of fourth order for a non-linear function f(R). Here prime (') denotes the derivatives with respect to R. 
The modified Eq. (fT3")) contains de-Sitter space time as a vacuum solution provided that /(4A) = 2A/'(4A). Thus, the 
f(R) theories of gravity may provide an alternative to dark energy. The f(R) gravity theories apart from a spin two 
object necessarily contain a scalar degree of freedom. Taking trace of Eq.fTjJ]) gives the evolution equation for the 
scalar degree of freedom, 

□/' = ~(2/'-/'#) + ^r. (w) 
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It would be convenient to define scalar function ip as, 

<f>=f-l, (15) 

which is expressed through Ricci scalar once f(R) is specified. We can write the trace equation (Eq. (fT4|) ) in the terms 
of V and T as 

(16) 



The potential can be evaluated using the following relation 

The functional form of f(R) should satisfy certain requirements for the consistency of the modified theory of gravity. 
The stability of f(R) theory would be ensured provided that, 

• f'(R) >0 — graviton is not ghost, 

• f"(R) > 0— scalaron is not tachyon. 

The f(R) models which satisfy the stability requirements can broadly be classified into categories: (i) Models in which 
f(R) diverges for R — > i?o where Rq finite or f(R) is non analytical function of the Ricci scalar. These models either 
can not be distinguishable from AC DM or are not viable cosmologically. (ii) Models with f(R) — > for R — > and 
reduce to cosmological constant in high curvature regime. These models reduce to ACDM in high redshift regime 
and give rise to cosmological constant in regions of high density and differ from the latter otherwise; in principal these 
models can be distinguished from cosmological constant. 

Models belonging to the second category were proposed by Hu-Sawicki and Starobinsky [r| [l7|. The functional 
form of f(R) in Starobinsky parametrization is given by, 



f(R) = R + XR 



1 + ' 



(18) 



Here n and A are positive. And i?o is of the order of presently observed cosmological constant, A = 8TrGp vac - The 
properties of this model [ItJ can be summarized as follows: 

1. Stability conditions are satisfied as /', /" > 

2. Flat space time is an unstable solution of the model. 

3. For \R\ ^> Rq, f(R) = R — 2A(oo). The high-curvature value of the effective cosmological constant is A(oo) = 
\Ro/2. 

In the Starobinsky model, the scalar field <p is given by 

4>{R) = m • (19) 

i?o(l + f)" +1 

Notice that R — > oo for <fi ~ * 0. One can easily compute V(R) for a given value of n. For instance, in case of n = 1, 
we have 

TT = — --^^4 {(-8-40y 2 -56y 4 -24y 6 )A+(32/ + lly 3 + 2l2; 5 -3y 7 )A 2 } 
-fto 24(1 + y 2 ) 

A 2 

-—tan" 1 y, (20) 
where y = R/Rq. In the FRW background, the trace equation (Eq. (fT4|) ) can be rewritten in the convenient form 

(j) + iH4> + — = — -p. (21) 
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FIG. 3: Plot of effective potential for n = 2 and A = 1.2. The red spot marks the initial condition for evolution. 



The time-time component of the equation of motion lfl~3|) gives the Hubble equation 



„ , d(ln/')„ , lf-f'R _ SnG 



H + v - ' H + -- — = p. (22) 

dt 6 /' 3/' F v ' 

It should be noted that the stability conditions ensure that the effective gravitational constant G e // — G/f appearing 
in Eq. (|22[) is positive. The simple picture of dynamics which appears here is the following: above infrared modification 
scale (Rq), the expansion rate is set by the matter density and once the local curvature falls below Rq the expansion 
rate gets effect of gravity modification. For pressure less dust, the effective potential has an extremum at, 

2/ - Rf = 8nGp. (23) 

For a viable late time cosmology, the field should be evolving near the minimum of the effective potential. The finite 
time singularity inherent in the class of models under consideration severely constraints dynamics of the field. 



The curvature singularity and fine tuning of parameters 



The effective potential has minimum which depends upon n and A. For generic values of the parameters, the 
minimum of the potential is close to <fi = corresponding to infinitely large curvature. Thus while the field is evolving 
towards minimum, it can easily oscillate to a singular point. However, depending upon the values of parameters, we 
can choose a finite range of initial conditions for which scalar field </> can evolve to the minimum of the potential 
without hitting the singularity. 

(p min 



-0.05 - 




-0.10 - 
-0.15 1 
-0.20 ; 
-0.25 1 



FIG. 4: Plot of (h m i n versus A for different values of n. With approaches zero for smaller values of A. The 

curves from bottom to top correspond to n = 1, 2, 3, 4 respectively. 

We find that the range of initial conditions allowed for the evolution of 4> to the minimum without hitting singularity 
shrinks as the numerical values of parameters n and A increase. This is related to the fact that for larger values of 
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n and A, the minimum fast moves towards <j> = 0, see figure [4] The numerical values should be accurately chosen to 
avoid hitting the singularity. 



Avoiding singularity with higher curvature corrections 

We know that in case of large curvature, the quantum effects become important leading to higher curvature cor- 
rections. Keeping this in mind, let us consider the modification of Starobinsky's model, 



f(R) = R+—R 2 + R Q \ 



-1 



1 



(24) 



then 6 becomes 



cj>{R) 



R 
Ro 



2a 



2n\ 



(1 



R 2 \n+1 



(25) 



When \R\ is large, the first term which comes from aR 2 dominates. In this case, the curvature singularity R = ±oo 
corresponding to <p = ±oo, see Figj5j Hence, in this modification, the minimum of the effective potential is separated 
from the curvature singularity by the infinite distance in the </>, V(<p) plane. In case of n = 2, the expressions for cj) 
and V(4>) are given by, 



</>(y) = 2ay 
V 

Rq 



AXy 



(l + y2)3' 



1 



480(1 + y 2 f 
1 



{X 2 y (-105 - 595y 2 - 2154?/ + 106?y 6 + 595z/ + 105y 10 )} 



. !( . + 8 - {1 + 5y 2 + a (3 + 8y 2 + 9y 4 + 4y 6 ) } + ±ay 2 + 1 (32a - 7A) tan^y). 



(26) 



(27) 



For n = 2, A = 2 and a = 0.5, we have a large range of the initial condition for which the scalar field evolves to 
the minimum of the potential. Though the introduction of R 2 term formally allows to avoid the singularity but 
can not alleviate the fine tuning problem as the minimum should be brought near to the origin to respect the solar 
constraints. Last but not the least one could go beyond the approximation (see Eq. (|2"5|) ) b y it erating the trace equation 
and computing the corrections to R given by equation (|23p . As pointed by Starobinskv [1 7| . such a correction might 
become large in the past. This may spoil the thermal history and thus needs to be fine tuned. The aforesaid discussion 
makes it clear that HSS models are indeed fine tuned and hence very delicate. 

In case of any large scale modification of gravity, one should worry about the local gravity constraints. The f(R) 
theories belong to the class of scalar tensor theories corresponding the Brans-Dicke parameter cu = or the PPN 
parameter 7 = (1 + lo)/{2 + uj) = 1/2 unlike GR where 7 = 1 consistent with observation (I7 — | < 10~ 4 ). This 
problem can be circumvented by invoking the so called chameleon mechanism. In case, the scalar degree of freedom 
is coupled to matter, the effective mass of the field depends upon the matter density which can allow to avoid the 
conflict with solar physics constraints. 



V. SUMMARY 



We have briefly summarized here various approaches to understand the late time acceleration of universe. In case 
of scalar field models of dark energy, we emphasized the relevance of scaling solutions in alleviating the fine tuning 
and the coincidence problems. We hope that the future data would reveal the metamorphosis of dark energy. 

Amongst the various alternatives to dark energy, the f(R) gravity models have received considerable attention in 
past years. There are broadly two classes of f(R) models, namely, those in which f(R) d iverges as R — ► Rq (Ro is 
finite) or f(R) is non-analytic in R. And those with f(R) — * as R — > 0, they reduce to KCDM in the limit of 
high redshift and give rise to cosmological constant in high density regime. These models can evade local gravity 
constraints with the help of the so called chameleon mechanism and have potential capability of being distinguished 
from ACDM^. 

Unfortunately, the f(R) models with chameleon mechanism are plagued with curvature singularity problem which 
may have important implications for relativistic stars (23|. The model could be remedied with the inclusion of higher 
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FIG. 5: Plot of the effective potential for n = 2, A = 2 and a = 1/2 in presence of R 2 correction, the minimum of the effective 
potential in this case is located at 4> m i n = 3.952 (Rmin = 3.958). 

curvature corrections [12] • At the onset, it seems that one needs to invoke fine tunings to address the problem f33[. 
The presence of curvature singularity certainly throws a new challenge to f(R) gravity models. In our opinion, the 
problem requires further investigations. It would also be interesting to look for a realistic scenario of quintessential 
inflation in the frame work of f(R) gravity. 
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